Semi- analytic Faddeev solution to the TV-boson problem with zero-range interactions 
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We study two-body correlations for N identical bosons by use of the hyperspherical adiabatic 
expansion method. We use the zero-range interaction and derive a transcendental equation deter- 
mining the key ingredient of the hyperradial potential. The necessary renormalization is for both 
repulsive and attractive interactions achieved with an effective range expansion of the two-body 
phase-shifts. Our solutions including correlations provide the properties of Bose-Einstein conden- 
sates exemplified by stability conditions as established by mean-field Gross-Pitaevskii calculations. 
The many-body Efimov states are unavoidable for large scattering lengths. 
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Introduction. Solutions to the two-body problem are 
described in many textbooks. The three-body problem is 
solvable in practice at least for short-range interactions 
PI Q- Analytic solutions are found for specific poten- 
tials, e.g. square- well potentials confined to only s-waves, 
where the large distance behavior is particularly simple 
. This limit is effectively the result for a zero-range in- 
teraction which is accurate when the small distances com- 
parable to the interaction range are uninteresting or only 
needed to provide overall average properties. Therefore 
much effort has been devoted to the three-body problem 
by direct use of zero-range interactions, where regular- 
ization is needed to prevent unphysical collapsed wave 
functions HHH- 

Not surprising the A-body problem has not been 
solved in general. The most popular approximation is 
the mean-field assumption which minimizes the techni- 
cal difficulties but ignores correlations. Even then of- 
ten the zero-range interaction is employed to exploit the 
additional simplifications allowing elaborate systematic 
investigations. Two prominent examples could be the 
Skyrme-Hartree-Fock fermion calculations for A-nucleon 
systems, see e.g.Q, and the analogous Gross-Pitaevskii 
calculations for dilute atomic Bose gases, see e.g. |3,|9(- 

Improvements within the mean-field approximation re- 
quire better interactions, for example of finite range. 
This changes significantly solutions of high density as 
for nuclei but would not be visible for dilute condensed 
atomic gases. Another improvement to include correla- 
tions and go beyond the mean-field is usually much more 
difficult. For nuclei a perturbative treatment starting 
with the zero-range interaction and the mean-field solu- 
tion is at the moment not very useful since a collapse 
into cluster components only can be avoided by a large 
additional phenomenological renormalization of the in- 
teraction. For dilute systems like atomic gases an ap- 
propriately renormalized zero-range interaction is not a 
severe approximation [lOj and it would therefore be di- 
rectly useful to include effects of correlations. 

Transparent or even analytic extension beyond the 
mean-field approximation is highly desirable for several 
reasons, e.g. (i) provide simplicity and insight, (ii) allow 
systematics and access of complicated systems, (iii) pro- 



vide intermediate and large distance properties, which 
in quantum mechanics often are directly responsible for 
the qualitative features of the entire solution, (iv) in de- 
signs of more efficient numerical methods by use of the 
available large distance asymptotic. The purpose of the 
present letter is to lay out the foundation for a series of 
applications and extensions by providing a semi-analytic 
solution to the TV-boson problem. 

Theory. We consider a dilute system of iV identical 
bosons of mass m, coordinates r, and momenta^. They 
affect each other via a two-body short-range interaction 
V(r) which can be approximated by a zero-range poten- 
tial. The Hamiltonian is then 
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where = \rj— r%\ and an external harmonic field of an- 
gular frequency uj is added. We choose the hyperspherical 
adiabatic expansion method where the principal coordi- 
nate is the hyperradius p defined by [l], [ll| 
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where R is the center of mass. The remaining degrees of 
freedom are the dimensionless hyperangles, collectively 
denoted by ft. The lowest adiabatic relative wave func- 
tion ^ is given by 
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where the angular part $ is expressed as a sum of two- 
body correlation amplitudes, Faddeev components, <f) 
each assumed to depend only on the overall size p and the 
distance between the pairs of particles. For zero-range 
interactions this dependence is no further restriction as 
two particles then only interact via s-waves. For fixed 
hyperradius the free variable in can conveniently be 
substituted by an angle atij defined by r = v^osinct, 
where we omitted the indices ij. 
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The Faddeev component (p(a) = 4>{p, vising) is de- 
termined from the angular Faddeev equation |llLll2j |. i.e. 



0= [ft 2 + v(a)R — X(p)] (j)(a) , 



(4) 



where h 2 X/(2mp 2 ) is the energy eigenvalue, v(a) is 
related to the two-body potential V(r) by v(a) = 
2mp 2 V(V2psina)/h 2 , and the kinetic-energy, ft 2 , and 
rotation, R, operators are given by 
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R = 1 + 2(N — 2)R 13 + i(JV - 2)(N - 3)R 34 . (6) 

The operators, Ri 3 and R34, in eq. © "rotate" two-body 
Faddeev components between particles 1 (or 2) and a 
third particle (2(N — 2) terms), and ones between parti- 
cles 3 and 4 both different from 1 and 2 ((JV-2)(iV-3)/2 
terms). The resulting radial equation is then 
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where E is the energy and the effective radial potential 
is 



2mU(p) _ A(p) (3iV-4)(3iV-6) p 
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with the oscillator length b 2 = h/(mu>). This potential 
is therefore determined by the external trap (p 2 ), the 
generalized centrifugal barrier (p~ 2 ) and the interaction 
part (A(p)) from the angular equation. The properties 
of the eigenvalue A from eq. (0} are then decisive for the 
radial potential and the corresponding energies and wave 
functions. 

Zero-range interaction. We approximate v by a zero- 
range interaction which means that eq. |0J has only ki- 
netic energy terms for all a / 0. The solution is well 
known as the Jacobi functions Vv (x), i.e. 

^(a) = ^ 3JV / 2 - 4 4/ 2 )(_cos2a) , (9) 
A = 2v{2v + 3N -5) , (10) 

where the boundary condition <j> u (a = tt/2) = is ful- 
filled by this Jacobi function HE)- boundary condi- 
tion of cj) v {a = 0) = is only obeyed for these functions 
for integer values of v which therefore fully determines 
the free solutions. However, at the point a = we now 
have an infinitely large potential which can be replaced 
by an appropriate boundary condition, i.e. obtained by 
using the observation that the wave function R<f) u (a) at 
small two-particle separation r approaches the two-body 
wave function u(r). Then v does not have to be inte- 
ger. The coordinate and wave function connections are 
r = V2p sin a w V2pa and u{r) oc aR(f) u (a). The bound- 
arycondition replacing the zero-range interaction for u 
is H 



1 du(r) 
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where k is the wave number, a s is the scattering length 
and R e ff is the effective range. For R e ff = we then 
get 
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For small a the solutions © behave as Q 
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Then at the edge of the zero-range potential we get 
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Combining eqs. (|12|l . (|15|l and (| 1 fc>|) we obtain 
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where (R — l)<f> v (0) are given by eq. © and the explicit 
expressions 
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The structure of eq. (|17fl is such that the right hand side 
depends on the particle number N and the index v on 
the Jacobi functions which in turn determines the angular 
eigenvalue A from eq.©. Thus the left hand side, p/a s , 
is a unique function of v and in turn of A through eq. IjlOj l . 
By inversion the effective radial potential is determined 
as a function of hyperradius divided by the scattering 
length. 

Extension to finite values of R e f / amounts to replac- 
ing l/a s in ea. H17f) by l/a s — ^k 2 R e ff, where mEijf? = 
k 2 = (A + (9JV - 19)/2)/(2p 2 ) with the two-body energy 
E2. This substitution is obtained by eq.(0J in the limit 
of small a from the connection between u(r) and </>„(«). 
Then eq. (|17f) becomes a second order equation in p where 
the physical solution easily is extracted. Now the effec- 
tive radial potential depends on both scattering length 
and effective range of the two-body interaction. 
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Angular eigenvalues. The expression in ea. (|17fl is 
straightforward to compute as a function of A via the 
(possibly complex) values of v. In fig. ^ we show the 
computed function p/a s as a function of A obtained as a 
sum of three contributions, i.e. the terms related to the 
Faddeev components between particles 1 — 2 (B-term), 
1-3 (eq. lfl%))) and 3-4 (eq. lO). This procedure is 
obviously easier than solving the transcendental equation 
to get A as a function of p/ a s . The immediate implication 
is that the interaction only enters via the ratio p/a s . 
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FIG. 1: The angular eigenvalues in units of Aoo = 
-1.65iV 7/3 (l - 2/N) as functions of p/(Na s ) obtained from 
eg, 1171 . The computation is for N = 100. The three different 
terms from eqs. I|17jl . i.e. B-term (dotted), 13-term in ea. l|18|l 
(long-dashed) 34-term in eq. 1191 (short-dashed), are shown 
individually along with the sum (solid). 

We confine ourselves to the most interesting cases of 
negative A corresponding to attractive two-body inter- 
actions. The terms 1 — 3 and 3 — 4 are both negative 
and vary from small A and large p/a s to large A and 
small p/a s . The first term is much smaller but crosses 
the p = axis and is therefore responsible for the same 
behavior of the sum of the three terms. However, this 
very important zero point for A is much larger for the 
sum than for the first term. We used units of Aqq from 
[l3l | and a scaling by N on the p-axis. Then the figure is 
fairly independent of particle number. 

When both the index and A are small, v <C 1, the Ja- 
cobi functions are almost constant allowing the analytic 
result 

^-i)^^ (20) 

where the validity condition, v « 1, can be translated 
into p ^> N 5 / 2 \a s \. This is the asymptotic behavior seen 
to the far left in fig. ^ F° r N ^> 1 where our center of 
mass separation is less important this result is identical 
to that derived in [ljj for a constant angular wave func- 
tion and a zero-range interaction renormalizcd for use in 
mean-field computations. 

At the other limit of large p/a s and large negative A 
we also obtain a closed analytic expression, i.e. 

Hp) - -2^J , (2i) 



where now only the first term in ea. (|17|l contributes. In 
fact this behavior corresponds precisely to the energy of a 
two-body bound state as E 2 = h 2 X/(2mp 2 ). We empha- 
size that the interaction only enters via the ratio p/a s . 
The scales in fig. ^ then implies that the results only ex- 
hibit the behavior for distances where p is (much) larger 
than a s . 
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FIG. 2: The angular eigenvalues in units of Aoo = 
-1.65N 7/3 (l - 2/N) for TV = 100 as functions of hyperradius 
for several values of the scattering length a 3 both in units of 
the effective range R e ff- 

For many realistic systems the interesting region is 
p/\a s \ < 1, i.e. the region on fig. ^ where p ps 0. To 
exhibit the behavior in this region we need to include the 
effective range in ea. <|llf) and solve the resulting second 
order equation in p/R e ff. The resulting angular eigen- 
values for different scattering lengths are shown in fig. [21 
as a function of hyperradius. No real solutions exist at 
small values of p where the detailed behavior of the two- 
body interaction in any case is important. Higher order 
terms in the effective range expansion would also allow 
the small p- values. 

As p increases the eigenvalue levels off when a s is suffi- 
ciently large. The height of this plateau is independent of 
a s and in fact precisely equal to the value of A obtained 
from fig.[2for P — 0. This value is rather accurately given 
by the unit used in fig.[T] A M = -1.65A^ 7 / 3 (1 - 2/N), ex- 
tracted by elaborate numerical solutions of a variational 
equation with finite-range interactions fl4| . 

As p increases further to values exceeding a s the eigen- 
values either bend up and approach zero (a s < 0) as ex- 
pressed in ea. (|2U|l or bend down diverging (a s > 0) as 
expressed in ea. l|21ll . These two characteristics then re- 
flect the different signs of the scattering length and the 
related ability of the two-body potential to support a 
bound state or not. 

Radial potential and solutions. The angular eigenval- 
ues are now inserted into the effective radial potential 
in eq.©- The behavior for different scattering lengths 
is seen in fig. [3] The confining trap provides the large 
positive potential for large distances. The minimum 
appearing for relatively weakly attractive potentials of 
\a s \/R e ff < 10.5 gradually disappears as the attraction 
increases. The barrier towards small distance disappears 
roughly when \a s \N/bt > 0.67 as derived previously in 



4 



0.003 



g 0.001 



-0.001 





~ \ ! < 1 

/ \ \ M / 
/ \ v 11 / 

/ " / 




/ ,»'"" a s /R eff =-2 

/ a s /R eff =-4 

:' , / a,/R. ff =-10.5 

/ / / a,/R B ff =-100 

J i non-interacting 



10 3 10 4 10 

0/R«£f 



FIG. 3: The effective radial potential in eq.@ for N = 100 
as a function of p for several values of the scattering length 
a s both in units of the effective range R e ff- The oscillator 
length is b% = 14427? e ff as in Hl| . The radial wave function 
(two-dashed line) located in the minimum is shown for a s = 
-2R eff . 

[l4l Il6| . The improvement is towards less stability de- 
creasing with increasing N. The minimum supports qua- 
sistationary states with characteristic features similar to 
condensates. The radial wave function / obtained from 
cq.Q) is shown in fig.[3]as a distribution around the min- 
imum. Other interaction parameters leaving the mini- 
mum would produce almost indistinguishable radial wave 
functions. 

The plateau region for large a s in fig. [21 now appears 
in fig. |3| as a p _2 -potential with a strength of approx- 
imately 2.25-/V 2 — 1.657V 7 / 3 , which is negative already 
when N > 3. As this strength is less than —0.25 the 
Efimov conditions are fulfilled and a number of states 
related by simple scaling properties are solutions to the 
radial equation in eq.JTJ), see 0, 0|. These states are 
located in the plateau region far outside the range of the 
two-body interaction but before the confining wall of the 



trap. If created they may have a sufficiently long life- 
time to be seen or perhaps play a role in some processes. 
They have obviously enough energy to decay into bound 
cluster states of much smaller size. 

Conclusions. Zero-range interactions have been used 
extensively for two and three-body systems. It is a sub- 
stantial simplification and very accurate for large dis- 
tances compared to the range of the potential. We extend 
the application to A^-boson systems by use of the adia- 
batic hyperspherical expansion method. We expect that 
two-body correlations are most important and dominated 
by s-waves. As zero-range interactions only are active in 
s-waves we use accordingly wave functions consisting of 
only s-wave Faddeev two-body amplitudes. This may 
then be viewed as the largest contribution in an expan- 
sion in both partial waves and many-body correlation 
amplitudes [jjj . For dilute systems we thereby obtain an 
accurate effective radial potential from a transcendental 
algebraic equation. 

We derive the crucial equation and renormalize the 
zero-range interaction by an effective range expansion in 
terms of two-body phase-shifts. This is analogous to the 
field theoretical renormalization studied intensively for 
three-body systems. We extract and discuss the perti- 
nent general scaling properties for both attractive and re- 
pulsive interactions. Use of parameters corresponding to 
systems forming condensates reveal the established prop- 
erties and stability conditions. 

The method has many applications and extensions 
which are beyond the scope of this letter, e.g. more 
systematic computations of various properties like one- 
and two-body densities, studies of correlations or effects 
beyond the mean-field approximation, investigations of 
dynamics in general and in particular the recombination 
process into bound cluster states, applications on more 
complicated systems like two-component boson systems, 
extensions to fermionic and mixed systems. 
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